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THE SIEGEL UPPER HALF SPACE IS A MARSDEN WEINSTEIN 
QUOTIENT: SYMPLECTIC REDUCTION AND GAUSSIAN WAVE PACKETS 

TOMOKI OHSAWA 


Abstract. We show that the Siegel upper half space is identified with the Marsden-Weinstein 
quotient obtained by symplectic reduction of the cotangent bundle with 0(2ci)-symmetry. 

The reduced symplectic form on Ed corresponding to the standard symplectic form on T E turns 
out to be a constant multiple of the symplectic form on Ed obtained by Siegel. Our motivation is to 
understand the geometry behind two different formulations of the Gaussian wave packet dynamics 
commonly used in semiclassical mechanics. Specifically, we show that the two formulations are 
related via the symplectic reduction. 


1. Introduction 

1.1. The Siegel Upper Half Space. The set of symmetric dxd complex matrices (symmetric 
in the real sense) with positive-definite imaginary parts, i.e., 

Srf := {A + iBe Mrf(C) \A'^ = A,B^ = B,B>0}, 

is called the Siegel upper half space. For d = 1, one easily sees that Si is the upper half plane 
El := {a -|- i5 G C | 6 > 0} of the complex plane C, and hence one may think of S^ as a generalization 
of the upper half plane H to higher dimensions. It is well known that the upper half plane El 
may be identified with the homogeneous space Sp(2,M)/U(l) via the linear fractional or Mobius 
transformation z i-G {az + b){cz + d)~^. The Siegel upper half space S^ is a natural generalization of 
El to higher dimensions in the sense that S^ is identified with the homogeneous space Sp(2d, M)/U(d) 
via a generalized linear fractional transformation; see Section [2] for details. 


1.2. Motivation: Gaussian Wave Packet Dynamics. Our motivation for studying the geom¬ 
etry of the Siegel upper half space is to better understand the underlying geometry for the 
dynamics of the Gaussian wave packet 


X{x) = exp 


i 

h 



q)'^{A + iB){x 


q)+p-{x 


q)+ {(!) + U) 


( 1 ) 


which is widely used in the study of the semiclassical limit of the Schrodinger equation. It is 
parametrized by {q^A G A + iB ^ f G and <5 G M, and it is well known (see Heller 

ji, and Hagedorn [J, [^) that ([I]) is an exact solution of the time-dependent Schrodinger equation 
with quadratic potentials if these parameters, as functions of the time, satisfy a certain set of ODEs. 

The geometry of —particularly the fact that is a symplectic manifold—becomes important 

when one tries to understand the set of ODEs as a Hamiltonian system on a symplectic manifold; 


see Ohsawa and Leok 


19| and Ohsawa 18l ]. 


Date: August 11, 2015. 

2010 Mathematics Subject Classification. 37J15, 53D20, 70G45, 81Q05, 81Q20, 81Q70, 81S10. 

Key words and phrases. Siegel upper half space, momentum maps, symplectic reduction, Hamiltonian dynamics, 
semiclassical mechanics, Gaussian wave packet. 


1 







2 


TOMOKI OHSAWA 


There are, however, two different ways of describing the dynamics. In the formulation originally 
due to Heller 0,3 j_^elements in are written as is, i.e., one writes down ODEs for A and B, 
whereas Hagedorn 0, 0] replaces ^ + iH by PQ~^ with d x d complex matrices Q and P that 
satisfy certain relationships, and the corresponding dynamics is written in terms of Q and P. The 
geometry behind the two different parametrizations turns out to be precisely the fact that 'Ed is the 
homogeneous space Sp(2d, M)/U((i), i.e., the set of variables {Q,P) naturally lives in the symplectic 
group Sp(2d, R) and ^ + iH is its projection to Sp(2d,M)/U(ci); see Sections [2] and 021 and also 
Ohsawa 


1 $. 


As simple as the correspondence sounds, one encounters an obstacle when trying to interpret 
the two formulations from the symplectic/Hamiltonian point of view. On one hand, it is fairly 
straightforward to formulate Heller’s dynamics with A and B from the symplectic/Hamiltonian 
point of view because is a symplectic manifold; see Ohsawa and Leok (l9l |. On the other hand, 
it is not so simple to do the same with the parameters Q and P of Hagedorn because the symplectic 
group Sp(2d,R) is clearly not a symplectic manifold in general; Its dimension is d{2d + 1), which 
is odd when d is odd. 


1.3. Main Results and Outline. The main result of the paper is Theorem 13.11 which is stated 
at the beginning of Section O In short, we show that the Siegel upper half space is identified 
with the Marsden-Weinstein quotient arising from the cotangent bundle with symmetry 

group 0{2d). Specifically, a certain level set of the momentum map and the corresponding isotropy 
group are identified with Sp(2(i,R) and U/d), respectively, thereby giving rise to the homogeneous 
space Sp(2d, R)/U(d) in the context of symplectic reduction. 

The theorem also gives a clear symplectic/Hamiltonian picture of the connection between the two 
different parametrizations of the Gaussian wave packet dynamics by showing that one is a Hamil¬ 
tonian system on T*R‘^ x with 0(2(i)-symmetry and the other is its reduced Hamiltonian 

system on r*R'^ x 

The outline of the paper is as follows: We first review, in Section [2l the geometry of the Siegel 
upper half space going through its realization as a homogeneous space; some of the results there 
will be later referred to when proving the main theorem. Section [3] states and proves the main 
result of this paper, Theorem 13.11 Finally, in Section01 we apply the theorem to the Gaussian wave 
packet dynamics to interpret the correspondence between two different formulations from the sym¬ 
plectic/Hamiltonian point of view. We also exploit the Hamiltonian nature of the problem to derive 
a semiclassical angular momentum by applying Noether’s theorem to semiclassical systems with 


rotational symmetry; this complements our earlier work on semiclassical angular momentum 181. 


2. The Siegel Upper Half Space Ed 

2.1. Geometry of the Siegel Upper Half Space. We first briefly review the well-known treat¬ 
ment of the Siegel upper half space Ed as a homogeneous space (see Siegel (20| and also Folland 
0, Section 4.5] and McDuff and Salamon [3, Exercise 2.28 on p. 48]). Specifically, we show the 
following identification alluded above: 


where Sp(2d, ] 


Frf^Sp(2d,R)/U(d), 

is the symplectic group of degree 2d over real numbers, i.e.. 


Sp{2d,R) := {S eM 2 d(M) \ S^JS = J} with J = 


0 

-Id 


Id 

0 
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or equivalently, written as block matrices consisting oi d x d snbmatrices, 


Sp(2d, 


'A B 
C D 


G M 


2d{ 


A^C = C^A,B^D = D^B,A^D-C^B = Id}, (2) 


and U(d) is the nnitary group of degree d. 

Consider the (left) action of Sp(2d, M) on defined by the generalized linear fractional trans¬ 
formation 


'k : Sp{2d, M) X Tid —)• '^d'l 
This action is transitive: By choosing 
X : = 


A B 
C D 


T” hA {C + DX){A + BA) 


-1 


A 

B 


'Id 

O' 

■i3-i/2 

0 


■ i3-i/2 

0 

C 

D 


.-4 

h 

0 

^1/2 


ab-^a 

i3i/2_ 


(3) 


which is easily shown to be symplectic, we have 

Tx(i/d) = ^ + iB. 

The isotropy subgroup of the element \Id G is given by 


Sp(2d,R)i/^ = 


U V 
-V u 


G M 


2d 


U^U + V^V = Id, u^v = v^u 


= Sp(2d,M) n0(2d), 


(4) 


where 0{2d) is the orthogonal gronp of degree 2d; however Sp(2d, M) n 0(2(i) is identified with U((i) 
as follows: 

■ U Fl 


Sp(2d,M) n0(2d) ^ U(d); 


-V U 


+ \v. 


(5) 


Hence Sp{2d,M)ii^ = U(d) and thns T^d = Sp(2d, R)/U((i). Indeed, we may identify Sp(2(i, R)/U(d) 
with Tid by the following map: 

Sp(2d,R)/U(d) ^ Brf; [T]u(d) ^ 'ky(i/d), 

where [ • ]u{d) denotes a left coset of U((i) in Sp(2(i, R); then this gives rise to the explicit construction 
of the qnotient map 


7ru(d): Sp(2(i,R) Sp{2d,R)/\J{d) = T ^'Y(i/d), 


( 6 ) 


or more explicitly, 


^U{d) 


A B 
C D 


= {C + [D){A + iB)-\ 


where A + iB can be shown to be invertible if £ Sp(2(i,R). 

As shown by Siegel [i^, B^ is eqnipped with the Hermitian metric 


:= tr {B-AC B-^dC) = Bji^B-\dCi^ ® dCjk 
and hence its imaginary part 

Hsd := ImfifSd = A dAjk = -dB~^ A dAjk 

defines a symplectic form on B^. 


( 7 ) 




















4 


TOMOKI OHSAWA 


2.2. The Iwasawa Decomposition and the Siegel Upper Half Space. We may make the 
above geometric structure more explicit by making use of the so-called Iwasawa decomposition of 
symplectic matrices (see, e.g., de Gosson (2, Section 2.2.2]). The Iwasawa decomposition renders 
any symplectic matrix Xi := £ Sp(2(i, M) into the following factorization of symplectic 

matrices: 


Xi = 


'A B 


1 

0 


'L 0 ■ 


■ U V 

C D 


P Id 


0 

1 


-V u 


L 0 
PL L-i 



■ u 

v 


-V 

u 


where P and L are symmetric d x d matrices given b’*^' 
p = {CA^ + DB^){AA^ + 
and X 2 and ZY G U (d) are defined b>l 


= X2U, 


L = {AA'^ + BB'^)^/'^, 


X2 : = 


L 0 
PL L-i ’ 


U:=U + IV = {AA^ + BB'^)-^/'^{A + iB). 


Then we see that the matrices Xi and X 2 define the same coset in Sp(2d,M)/U(d), i.e., [Xi]u(d) = 
[-^2]u(d) because o ^uiih) = ^X 2 ^'^h)- 

Now, take any element .4. -|- i.S in the Siegel upper half space and suppose that 7ru(rf)(Xi) = 
A + iB. Then 7ru((i)(Xi) = 7ru(d)(X2) = P + iL~^ using ([3]) and so P = A and L = B~^l‘^. 
This observation leads to the following expression for those elements in Sp(2d, M) that project to 
A + iB £T.d via 7ru(d): 




■ i3-i/2 

AB-^/^ 


0 

51/2 



■ u 

v 


-V 

u 


G Sp(2d,M) I U + iU G U(d) 


( 8 ) 


This expression will be later useful in coordinate calculations. 


3. The Siegel Upper Half Space and Symplectic Reduction 


In this section, we prove our main result, that the Siegel upper half space is identified with 
a reduced symplectic manifold in the Marsden-Weinstein 1^ sense. More specifically: 

Theorem 3.1. Let := he the cotangent bundle of the vector space with the standard 

symplectic form where we identify with the set Mrfx 2 d(I^) of d x 2d real matrices and 

Zd with the set M 2 (i(M) of 2d x 2d real matrices. Consider the action of the orthogonal group 
0{2d) on Zd defined by matrix multiplication from the right, and let M: Zd ^ o{2d)* be the 
corresponding momentum map. Then the Marsden-Weinstein quotient Zd '■= lVI~^(J)/0(2d)j at 


J = 


0 Id 
-Id 0 


G o(2d)* is the Siegel upper half space with symplectic form Llj = — i.e. 


where ij: Zd is the inclusion and ttj: — )• Zd is the projection to the quotient. 


We prove the above theorem in the rest of the section. 

^The matrix AW + BB^ is always invertible if Xi is in Sp(2ci,R) and hence so is L; see, e.g., de Gosson [3, 
Section 2.1.2]. 

^Note that, in the above expression X 2 U, the matrix U is seen as an element in Sp(2d,R) n 0(2d) via (]5|l. 
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3.1. Basic Setup. Consider the real vector space we prefer to write each element in as 
a d X 2d block matrix consisting of two d x d submatrices, i.e., 

= Mrfx2d(M) := { [Qi Q 2 ] I Qi,Q 2 G Md(M)}. 

Let Zd := be the cotangent bundle of Then each element Z in Z^ is identified with a 

2d X 2d block matrix consisting of four d x d submatrices, i.e., 

Zd := ^ M2d(M) = |z := 

The standard symplectic structure on Zd is given by 

Qza •= dQi A dPi + dQ2 A dP2 = dQ{^ A dPijk + dQ^^ A dP2,jfc, (10) 

where j and k run from 1 to d, and and Pijk stand for the (j, A;)-entries of the matrices Qi and 
Pi, respectively, for / = 1,2. With the canonical one-form Qz^ on Zd defined by 

Qza ■= tr(pfdQi) + tr(P 2 ^dQ 2 ) = Pi,jkdQ{^ + P 2 ,jkdQf, 
one can write the symplectic form ^Zd 


Qi Q2 
Pi P2 


Qi,Q2,Pi,P2 G Mrf(M) 


(9) 


^Zd = -^P>Zd- 


3.2. O(2(i)-action and Momentum Map. Consider the action of the orthogonal group 0(2d) on 
Zd defined by matrix multiplication from the right, i.e.. 


^■.0{2d)xZd^Zd-, {n,z)^zn. (ii) 

It is a straightforward calculation to see that $ leaves the canonical one-form Qzd invariant, i.e., 
^Tl^Zd = ®Zd for any TZ € 0{2d), and hence is a symplectic action with respect to the symplectic 
form ^Zd defined in (flUjl . i.e., = Uz^ for any TZ E 0(2d). 

What is the momentum map corresponding to the O(2(i)-action <I>? Let o(2d) be the Lie algebra 
of 0(2d) and o{2d)* be the dual of o{2d). For any ^ E o{2d), the corresponding infinitesimal 
generator ^Zd is given by 


^ZdiZ) := ^^exp(£5)(^) 




£ = 0 


where Z^ stands for the standard matrix multiplication. Then, since Zd is an exact symplectic 
manifold with the symplectic structure Qzd = —dOzd the action leaves 0Zd invariant, the 
corresponding momentum map M: Zd —>■ o(2d)* may be defined as follows (see, e.g., Abraham and 
Marsden [l|. Theorem 4.2.10 on p. 282]): For any ^ E o(2d), 


(M(Z),0 = Ozd(^Zd(Z)) or (M(-),0=k^^eZd. 
We equip the Lie algebra o(2d) with the inner product 


(•, •): o(2d) X o(2d) 


{^,V) ^ (Lv) ■= ^tr(C'^r/). 


( 12 ) 


Then we may identify the dual o(2d)* with o{2d) itself via the inner product. So we may write the 
components of the momentum map M as follows: 


M: Zd^o{2d)*] 


Z = 


Qi Q2 
Pi P 2 


I-;- 


Mil Mi2 

—MJ2 M22 


where both Mu and M 22 are skew-symmetric d x d matrices and M 12 is a d x d matrix, i.e., 
Mu,M 22 E o(d) and M 12 E Md(M). 
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Let US first find Mu. Pick ^ G o{2d) with G o(d); then 


(M(Z),e) = ^tr 


^11 0 
MJ 2 C 11 0 


1 


= -tr (Mfi6i). 


Qi 

Q 2 

All 

0 ’ 


Qi Cii 0 

Pi 

P 2 

0 

0 


Pi 61 0 


On the other hand, the corresponding infinitesimal generator is given by 

^zAz) = 

and hence 

ezACzAz)) = tT{PAQiCii) = ^tr [(Q^Pi -pfQifeii]. 
Since G o{d) is arbitrary, we find Mu = QjPi — P^Qi- 


Likewise, ^ 

Finally, taking ^ = 


G o{2d) with ^22 G o{d) yields M 22 = Q 2 P 2 - P^Q^- 
G o{2d) with I 12 G we have 


0 ?12 
-sr2 0 


1 


{MiZ),0 = ^tT 


Mi2 ^^2 
-^^2 ^12 


mJ[ ^12 
MT 2 A 2 


= tr {MJ 2 62 ), 


Qi 

Q 2 

0 

62 


Q 2 ?12 

Qi 62 

Pi 

P 2 

1 - 

1 

0 


rP^iT2 

Pi ii2 


whereas the corresponding infinitesimal generator is 

^zAz) = 


and so 

ezA^zAz)) = tii-p[Q 2 ^l 2 ) + tAPAQiA 2 ) = tr [(gfP 2 - PAQ 2 fA 2 ]. 
Again, since ^12 G Mrf(]R) is arbitrary, we find M 12 = Q 1 P 2 — Pi^g 2 - 
As a result, we have the momentum map 


M: Zw 


o(2d)*; 


Z = 


Qi Q 2 
Pi P 2 


i-A 


QJPi-PAQi 

-iQjP,-pTQ,)T 


QIP 2 - pTQ2 
QIP 2 -PAQ 2 


(13) 


It is a straightforward calculation to check that M is equivariant, i.e.. 

Mo $ 7 ^ = Ad^oM. 


3.3. Momentum Level Set and Reduced Space. Now let us look at the level set of the mo¬ 


mentum map M at J = 


0 Id 
-Id 0 


G o{2d)*, i.e.. 


M-i(J) 


Qi 

Pi 


Q 2 

P 2 


G M2d(M) I QfPi = PfQi, QIP 2 


PAQ2, QIP2 - PIQ 2 = Id 


This is precisely the definition of the symplectic group Sp(2(i, M) in terms of block matrices shown 
in ([2]), i.e., M~^(J) = Sp(2(i,M). 

Moreover, the coadjoint isotropy subgroup 0{2d)j of J G o(2d)* is easily identified as 

0(2d) j = {P G 0{2d) I Ad^ J = PAJU = J} 

= Sp(2d,M) n0(2d) ^ U(d), 


which is precisely the isotropy subgroup Sp(2d, R)i 7 ^ of the action of Sp(2(i, M) on the Siegel upper 
half space shown in Q. The action of the coadjoint isotropy subgroup 0{2d)j on the momentum 
level set = Sp(2d, M) is free, as STZ = 5* for 5 G Sp(2(i, M) and TZ G 0{2d) implies IZ = l 2 d'i 

the action is also proper as well because U((i) is compact. 
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So we may now invoke the Marsden-Weinstein reduction 0 (see also Meyer 0] and Marsden 
et al. [la, Sections 1.1 &: 1.2]) to conclude that the reduced space or the Marsden-Weinstein quotient 
:= M“^(J)/0(2d)j is a symplectic manifold, but then this quotient coincides with the Siegel 
upper half space = Sp(2d, M)/U(d), i.e., 

Zi := M-i(J)/0(2d)j = Sp(2d,M)/U(d) ^ 


3.4. Reduced Symplectic Form. Let us define the following inclusion and projection maps: 


ij: M-\J) -A Zrf, vrj: M-'(J) ^ M-'(J)/0(2d)j =: 

As shown in Marsden and Weinstein 1^, the symplectic form 11 j on the reduced symplectic 
manifold Z^ corresponding to the original one is uniquely characterized as follows; 


7T*jQj = i*jnz,. (14) 

We would like to find an expression for 11 j. Let ^ + i.8 be an arbitrary element in Since 

7ru(rf) in ([6]) and ttj are identical, those elements in Sp(2(l, M) that project to A + iB are written as 
in ([8]). But then this implies that any element in Sp(2d,M) = 1V[“^(J) = 7rJ^(E^) is written as 


■ ^-V 2 

0 

' u v' 



^-i/ 2 y 

AB-^/^ 

^ 1/2 

-V u 


AB-^/^U - B^/‘^V 

AB-^I'^V + 


with some A + iB E and U + iV £ U((l), i.e., U'^U + V'^V = Id and U'^V = V'^U. In 
other words, the above expression gives an expression for the inclusion ij: M“^(J) Zd in 
terms of the coordinates adapted to the horizontal and vertical directions of the principal bundle 
TTj: M“^(J) —)• M“^(J)/U(d). Then the pull-back by ij of the one-form Qz^ is written as 

i*jQz^ = ]^ti{AdB-^) + ti{U^dV -V^dU), 

and so, taking into account the relationships between U and V, 

i*jllz, = -d(i}02j = ^dB-,1 A dAjk, 

and hence we have, from (fT4]) and dT]), 7r}llj = — Since ttj is a surjective submersion, tt} 

is injective; thus we obtain Hj = — This completes the proof of Theorem 13.11 


4. Application to Gaussian Wave Packet Dynamics 


4.1. The Gaussian Wave Packet Dynamics. Consider the time-dependent Schrodinger equa¬ 
tion 


d 

ih—MxA) = AV’(x,l) + V{x)i;{xA) 
ot 2m 


(15) 


for the wave function 'il){x,t) under the potential P(x), where h > 0 is Planck’s constant, 1 > 0 is 
the time, x E is the position in the physical space M'’*, and A stands for the Laplacian in 
Our motivation for identifying the Siegel upper half space as a Marsden-Weinstein quotient 
comes from a geometric description of the dynamics of the Gaussian wave packet ansatz 

/detf i r 1 


V’o := 


V(7rh)'^y 


exp 


(x - qY {A + iB){x - q) + p ■ {x - q) + (p 


(16) 


for (|15p : the factor in front of the exponential normalizes the wave function, i.e., ||V’o|| = 1 as an 
element in L^(]R‘’*), and Yo is parametrized by {q,p) E C '■= A + iB £ T^d, and (p £ S^. It is 
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well known (see Heller [g, L^) that, when V is quadratic, the Gaussian wave packet (fT6t) is an exact 
solution of the Schrodinger equation (fTSl) if the parameters {q,p,A,B) satisfy the set of ODEs 

p_ 

m ’ 

22\ V72t 


9 = 


P = -VH(g), 


1 


1 


(17) 


A = - {A^ -B^)-V^V{q), B = - (AB + BA), 


m 


m 


where is the Hessian matrix of V, and the phase (/>(t) is determined by 


= 4>{o) + [ 

Jo . 




ds. 


Hagedorn [J, l5||, on the other hand, has a slightly different parametrization of the Gaussian wave 
packet (fTBI) : 

1 


= (7rh)-'^/^(detQ)-^/2 


exp 


h 


{x - qy PQ yx - q) + p ■ {x - q) + S 


(18) 

(19) 


where Q and P are complex dx d matrices, i.e., Q,P G Md(C), that satisfy 

qTp _ pTq ^ 0 and Q*P - P*Q = 2ild, 
and an appropriate branch cut is taken for (det also the new parameter S is defined as 

S := (p- ^arg(detQ). 

Hagedorn [3, showed that (fTS|l is an exact solution of the Schrodinger equation if the potential 
V is quadratic and also the parameters {q,p,Q,P) satisfy 

P 


Q = 


( 20 ) 


m m 

and the quantity S{t) is the classical action integral evaluated along the solution (g(t),p(t)), i.e., 


S{t) = S{0) + l^ (^P^-V{q{s))ys. 


It is also shown by Hagedorn [J, that (fTSl) with (1201) gives an 0{ty/h) approximation when the 
potential V is not quadratic as long as it satisfies some regularity assumptions. 

4.2. Parametrization of the Siegel Upper Half Space. The replacement oi A + \B G T^d in 
(fT6l) by PQ~^ in (fTSl) has a simple geometric interpretation. Lubich [l^, Section V.l] (see also 
Ohsawa [ist]) pointed out that the conditions (IlOp for the matrices Q and P are precisely the 


conditions for the matrix 

Sp(2(i, M) = 


Re Q Im Q 
Re P Im P 


Re Q Im Q 
Re P Im P 


In fact, the projection of these elements to by : Sp(2(i, 

Re Q Im Q 


to be symplectic, i.e., 

Q,Pg Mrf(C), Q^P-P^Q = 0, Q*P-P*Q = 2i/rf|. 

Erf in dH) gives 


^U(d) 


Re P Im P 


= PQ 


-1 


( 21 ) 


It is also easy to show that the dynamics of A and B defined in (fT7)l is the projection to by 7ru(d) 
of the dynamics of Q and P in (p0|) : conversely, the dynamics ([20]) is a proper lift to Sp(2(i,M) of 
the dynamics (fT7)l in some appropriate sense; see Ohsawa [l8|. 
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4.3. Interpretation as a Hamiltonian Reduction. Theorem 13.11 sheds a new perspective on 
the above connection between the equations (1171) and (|20p of Heller and Hagedorn, respectively, in 
terms of Hamiltonian reduction. 

Let us first introduce a new parametrization of the space Z^: The observation (|21l) from the 
previous subsection motivates us to rewrite the definition Q of the space Z^ as follows: 


:= T* 




= lz = 


Re Q Im Q 
Re P Im P 


Q,Pe Md(C) 


i.e., Qi + \Q 2 = Q and Pi + iP 2 = P. 

Here we assume that the potential V is quadratic for simplicity, and will consider the general 
case in the next subsection. If V is quadratic then the Hessian is a constant matrix; as a 

result, the system (flT]) decouples into the classical Hamiltonian system in r*M'^ = {{q,p)} and 

A = - — {A^ 13 = - — {AB + BA) ( 22 ) 

m m 

in Zd = {(^,13)}, and similarly, (I2nh decouples into the classical Hamiltonian system in = 

{{q,p)} and 

Q = -, P = -V^VQ (23) 

m 

in Zd = {{Q, P)}. It turns out that the dynamics (l2^ in is obtained by Hamiltonian reduction 
of the dynamics (l23l) in Zd‘. 


Proposition 4.1 (Reduction of Gaussian wave packet dynamics—quadratic potentials). Suppose 
that the potential V is quadratic, and let Hz,i: Zd ^ ^ be the Hamiltonian defined by 

HzfiZ) := ^ tr(P*P) + i tr(g* VV Q). (24) 

Then: 

(i) The Hamiltonian vector field on Zd defined by the Hamiltonian system 

ixn (25) 

gives the equations (l23l) of Hagedorn. 

(a) The Hamiltonian Hz^ is invariant under the 0(2d)-action $ defined in (fTT]l . i.e., Hz^ o 
= Hza for any TZ € 0(2d), and hence the Hamiltonian system (|25]) conserves the 
corresponding momentum map M in (HSI) ; particularly, M ^(J) = Sp(2(i, M) is an invariant 
manifold of the system (l25l) . 

(Hi) The reduced Hamiltonian Hj\ Z^^ M, i.e., the function Hj uniquely characterized by 


takes the form 


HjoTTj = Hza oij, 


Hj{A,B) = 2 




m 


(26) 

(27) 


(iv) The Hamiltonian vector field on the reduced space Z^d — defined by the Hamiltonian 
system 

ixjj = dHj 

gives the equations (1221) of Heller. 


(28) 
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Proof. The first assertion follows from straightforward calculations; so is the 0(2(i)-invariance of 
the Hamiltonian ()24p . 

The remaining assertions follow easily from Theorem 13.11 By Noether’s Theorem (see, e.g., 
Marsden and Ratiu IJ, Theorem 11.4.1 on p. 372]), the 0(2(i)-invariance of the Hamiltonian Hz^ 


implies that the momentum map M is conserved along the flow defined by Xhz ; we note that this 


I— 


result is observed by Lubich 113 , Lemma V.1.4 on p. 126] via direct calculations. We have already 
shown in Section [T3] that M“^(J) = Sp(2d, M). 

The reduced Hamiltonian (12711 follows from the defining relation (I26p and coordinate calculations 
that are similar to those performed in Section [3.41 when finding the reduced symplectic form Hj. 
Finally, that the reduced Hamiltonian system is given by (|28l) is standard in symplectic reduc- 

□ 


tion 


141 ]. and straightforward calculations yield the last assertion. 


Remark 4.2. For the Hamiltonian dynamics ()25ll to be interpreted as the dynamics of the Gaussian 
wave packet (|18ll . one needs to restrict the initial condition to the invariant manifold M“^(J) = 
Sp(2d, M) C Zd, which is equivalent to (|19p as mentioned above. This guarantees that Q is invertible 
and also that PQ~^ G and hence fjo G 


Remark 4.3. Setting C := A + iB, (I22p is written as a matrix Riccati equation, i.e., 

c = - vV. 

m 

Its Hamiltonian lift (|23p to Sp(2d, M) C Z^ is linear in Q and P, and may be considered as an 
example of the Hirota bilinearization of the matrix Riccati equation; see, e.g., Hirota 


4.4. Pull Dynamics as Hamiltonian Systems. What if the potential V is not quadratic? It 
turns out that, with a slight modification in the equations (|17p and (|20p . these systems may also 
be rendered Hamiltonian as well, and again one is the reduced version of the other. 

Let ■— ^9* dpi be the standard symplectic form on pr^: x —)• T*M.^ and 

pr 2 : X Zd ^ Zd be the projections, and define a symplectic form fl on T* X Zd by 

H := pr]' + — pr^ 11 (29) 

and a Hamiltonian P: r*]R'^ x —)• M by 

= L- + T tr(P'P) + V(q) + 1 tr((3" VV Q). (30) 

2m 4m 4 

Note that the symplectic form 17 is written as 17 = —d0 with 

0 := pi’i ©r.ffid + ^ pr 2 Qza = Pj ^ [tr(Pf dQi) + tr(P 2 ^dQ 2 )] • (31) 

Proposition 4.4 (Reduction of Gaussian wave packet dynamics — general potentials). Let V G 
C'^(R'’*) and H: x —>■ M be the Hamiltonian defined in (|30ll . 

(i) The Hamiltonian system 


ix^H = dH 


(32a) 
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gives the system 

# = -^{u9) +jtr|(3*vT(,)<3]}. 

Q = -, P^-VH'(g)Q. 

m 


(32b) 


(ii) The Hamiltonian dSO]) is 0{2d)-invariant under the 0{2d)-action on r*M'^ x 

and hence the corresponding momentum map M: r*M'^ x Zd ^ o{2d)*, which is given by 
M = M o pi 2 , is conserved along the flow of the system ([32]). In particular, = 

r*M'^ X Sp(2d,M) is an invariant manifold of (1321) . 

(Hi) Symplectic reduction by the 0{2d)-symmetry at the value J E o{2d)* yields the reduced 
symplectic manifold T*]R'^ x with symplectic form 


h 


It := dg* A dpi + - dBjif A dAjk 


and the reduced Hamiltonian 
12 


H = ^ + V{q) + ^ir 
2m 4 


B~^( ^^-^ + VV(g) 


m 


(33) 


(34) 


i.e., they are uniquely determined byXfjH = IIjO and HoXj = HoHj, whereXj: M ^(J) 
r*M'^ X Zd and IIj: M“^(J) —)• T* X Zd are the inclusion and projection, respectively, 
defined by 

Xj '■= id'jitjjd xij, IIj := idj’tRd xttj. 


(iv) The reduced Hamiltonian system 


ix-It = dH 


(35a) 


gives 


P 

q = 

m 
1 


d 


V{q) + -tT (^-ivV(g)) 


(35b) 


A = --{A^-B^)-V^V{q), 
m 


B = - (AB + BA). 

m 


(v) Particularly, when the potential V is quadratic, the systems (|32|) and (135|) recover (120|) and 
m, respectively. 

Proof. Most of the assertions are straightforward generalizations of those in Proposition 14. II and so 
are proved in a similar way. The last assertion follows because the Hessian becomes constant 
when the potential V is quadratic. □ 


Remark 4.5. The above symplectic form Q and Hamiltonian H are exactly those that appeared 
in Ohsawa and Leok [I^ (see 11^ in (17) of Theorem 4.1 and also Hi in Section 7.1 there; see 
also Ohsawa 1^, Section H.C]), but are derived from a different point of view. It was shown in 
Ohsawa and Leok 1^ that the symplectic form 0 in (I33p is a pull-back of the symplectic form on 
the projective Hilbert space P(L^(M'^)) induced by the symplectic form 


^^(V’U'02) = 2;iIm(V’i,V'2) 
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on —the Schrodinger equation (11511 is written as a Hamiltonian system on or P(L^(M'^)) 

in terms of those symplectic forms; see, e.g., Marsden and Ratiu Ij, Chapter 2]. As for the Hamil¬ 
tonian H in dMl), with a certain decay property assumed for the potential V, one can show that H 
is an 0{h^) approximation to the expectation value of the Hamiltonian operator 


H:=-^A + Vix) 


with respect to the Gaussian (fT6]l . i.e., 




where (•, •) stands for the inner product on see Ohsawa and Leok 19|, Proposition 7.1]. 


4.5. Rotational Symmetry and Conservation Law. The 0{2d) symmetry exploited above is 
rather an intrinsic symmetry of the Hamiltonian system ()32p . Here we would like to see what hap¬ 
pens if the system has a symmetry in its physical configuration, particularly when the potential V 
has a rotational symmetry. This leads to the semiclassical angular momentum for the semiclassical 
system (f32]) . which corresponds to the one for the reduced system ([3^ found in Ohsawa [l8| : 


Proposition 4.6 (Rotational symmetry and semiclassical angular momentum). Let (p: SO{d) x 
be the natural action of the rotation group SO{d) on the configuration space R'^, i.e., for 
any R E SO{d), 

ipn: ^ qe^Rq, 

and suppose that the potential V E C'^(R'^) is invariant under the SO{d)-action, i.e., 

V o PR = V. 

Also let T: SO{d) x (T*W^ x Z^) —>• T*W^ x Z^ be the SO{d)-action on T*W^ x Z^ defined as follows: 
For any R E SO{d), 

Tr : T^R-" X Zrf ^ r^R'^ x Zp, {q,p, Q, P) ^ {Rq, Rp, RQ, RP). 

Then: 

(i) T leaves the canonical one-form © invariant, and hence is a symplectic action with respect 
to the symplectic form (12^ . i.e., TrLI = for any R E SO{d). 

(a) The Hamiltonian ([30]l is invariant under the action, i.e., HoTr = H for any R E SO(d). 
(Hi) The semiclassical system dMl) conserves the semiclassical angular momentum J: T*R'^ x 
Zd —)■ 5o{d)* defined by 

3{q,p,Q,P) =qop+^{PiQl + P^QI - QiPf - Q 2 PI), (36) 

where qop denotes the dx d matrix defined by {qop)ij ■= qjPi — qtpj; see, e.g.. Holm fjh. 
Remark 6.3.3 on p. 150]. The angular momentum map J is equivariant as well, i.e., for 
any R E SO{d), 

3oTr = J. 

Proof. R is straightforward computations to see, for any R E SO{d), that T^0 = 0 and so 
HrLI = 0 as well as that H o Tr = H. For the symplecticity of T, one may alternatively identify 
T*R'^ X Z^ with the cotangent bundle r*(R'^ x M(^x 2 (i(I^)) and then see that T is the cotangent lift 
of the action 

SO{d) X (R*^ X Mdx 2 d) ^ 0^'^ X Mrfx 2 d); {R, (q, <3i, Q 2 )) {Rq, RQi, RQ 2 ); 
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see, e.g., Marsden and Ratiu [ij, Proposition 6.3.2 on p. 170]. 

Let us find the corresponding momentum map, i.e., the angular momentum for the semiclassical 
system (I32|] . For any ^ € so(d), its infinitesimal generator on x is given by 

d 




d 


d 


e=0 

2 


d 


= +E ( 


d 


dp 


1=1 


dQ 


dP, 


l,jk 


where the indices j and k run from 1 to d. Since T*M'^ x Zj^ is an exact symplectic manifold 
with symplectic form fl = —d0 with 0 given in ()31l) . and T leaves 0 invariant, the corresponding 
angular momentum J: x Z^ ^ 5o{d)* is given by 


(J(g,p,Q,P),0 = or (J(-),0 = 

and straightforward computations yield 

{ 3 {q,p,Q,P ),0 =P-iq + |tr(Pf^Qi + P2 iQ2) 

= (qop+^{PiQi+P2Ql - QiPf - Q2P^),^y 

where we identified 5o{d)* with 5o{d) via the inner product as in (|12p : 

(•, •): 50{d) X 50{d) M; {p rf) ^ {p 77) := ^ tiif'r]). 

So we obtain ([36]1 : it is a conserved quantity of the system (1321) due to the SO((i)-invariance of 
the Hamiltonian H and Noether’s Theorem (see, e.g., Marsden and Ratiu [lil . Theorem 11.4.1 on 
p. 372]). 

The equivariance of J is easy to show by direct calculations: 

J{Rq, Rp, RQ, RP) = R 3 {q,p, Q, P)R^. 


The equivariance also follows from the fact that T is a cotangent lift as mentioned earlier; see, e.g., 
Marsden and Ratiu [3, Theorem 12.1.4 on p. 386]. □ 


Assuming that {q,p,Q,P) is in the invariant manifold M ^(J) = r*M“xSp(2d, M), i.e.. 


Re Q Im Q 
Re P Im P 

Sp(2(i,M) or equivalently (fl^ . the setup and result in the above proposition descend from 
to X and recover the semiclassical angular momentum found in Ohsawa 18]. In fact, the 

action T induces an SO((i)-action T: SO{d) x (T*W^ x —)• T*M'^ x so that the diagram 


T*R'^ X Sp(2d,l 
nj 

X Erf - 

commutes: For any R G SO{d), we have 
Fr : r*R^ X Erf r*R‘^ x Eq 


r*R'^ X Sp(2fi,] 


A T* 


nj 

X Ef^ 


{q,p,A,B) I—)• {Rq, Rp, RAR^, RBR^). 


(37) 


This coincides with (26) in 0hsawa 1^; in fact F is a symplectic action with respect to the 
symplectic form (l33]l . i.e., FJ^H = Q for any R € SO(d). 
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Notice that 

Re(Pg* - QP*) = PiQl + P 2 Q 2 - QiPl - Q 2 PI 
with Q = Qi + iQ 2 and P = Pi + \P 2 as defined above. Then (fT^ and (1211) give (recall that 
'^U(d) = T^j) 

ttj{Q,P) = PQ-^ = A + 'iB 

as well as QQ* = because (Lubich [H, Lemma V.1.1 on p. 124]) 

BQQ* = Im(Pg-i)QQ* 

= -{Q*)-^P1QQ* 

= ^^{PQ* - {Q*)-^P*QQ*) 

= ^ W - {Q*)-\Q*P - 2ih)Q*] = h, 

where we used the second equality of (fT^ . Therefore, 

PQ* - QP* = PQ-\QQ*) - Q{P*Q)Q-^ 

= PQ-\QQ*) - Q{Q*P - 2i/rf)Q-i 
= (A + iB)B-^ - B~\A + iB) + 2ild 
= [A,B-^] + 2iId, 


and hence we have J 


r*R'*xSp(2d,R) 


= J o IIj with J : r*M‘^ x 5o{d)* defined by 


J{q,p,A,B) = qop-^[B ^,A]. 


This is exactly the semiclassical angular momentum for (I35|) derived via the action (I37p in Ohsawa 


la . Theorem 3.1]. 


The above semiclassical angular momentum is a natural one in the quantum mechanical sense 
as well, as mentioned in Ohsawa 18], Section III.B]. In fact, set d = 3 and let x and p := —i/iV be 
the position and angular operators; then one can show that, again assuming {q,p, Q, P) E 

J{q,p,Q,P) = '^{q,p,AP) = (V’o, {x X p)i^o), 

where ^po is either (|16p or m- 
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